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Abstract 

The phenomenon of quantum tunneling is reviewed and an overview of applying approx- 
imate methods for studying this effect is given. An approach to a time-dependent formalism 
is proposed in one dimension and generalized to higher dimensions. Some physical examples 
involving the resulting wavefunction which is determined are presented. 
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There is a great deal of interest in integrable systems at the moment and their quantization. 
The semiclassical treatment of integrable systems has its roots in Bohr's atomic model as well 
as in Einstein's well known paper on the quantization of regular motion [1]. It is known that 
semiclassical quantization is exact only for so-called solvable potentials, which would include 
the harmonic oscillator [2,3]. The nonsolvable case yields approximate results which may be 
asymptotic in nature [4,5]. Of course there are summation techniques which can be applied 
to divergent series. Recently [6] a quantum phase was derived which gives exact quantization 
of nonsolvable potentials and avoids many of these complications. A consistent divergence-free 
scheme is still an active area of research. 

The phenomenon of tunneling has been a characteristic property of many types of quantum 
mechanical processes, especially those which involve confining potential wells [7]. Their math- 
ematical description has generated many different types of results. Due to the wide variety of 
applied applications where tunneling plays some role, the study of mathematical methods which 
would be useful in this context is still very active. Along these lines, the quantum WKB approxi- 
mation is very important, and not a great deal exists in the literature for the time- dependent case 
[8]. An important example of a system where this approximation can be useful is the case of the 
decay of a system from a long-lived metastable state which frequently arises in physics. There 
has always been considerable interest in understanding the quantum tunneling of a macroscopic 
coordinate. If the height of a potential barrier is greater than the total energy E of a particle, the 
kinetic energy of the particle is negative and so the study of the tunneling of a particle appears to 
have a paradoxical aspect to it. However this approach implies that we can assign values to the 
coordinate x and the momentum p simultaneously [7]. This however would violate the uncertainty 
principle. If however for a very short time At, the uncertainty in the energy AE is such that the 
total energy of the particle is greater than the height of the barrier, then tunneling takes place 
in time At if the particle can traverse the barrier in this time interval. The continued interest in 
the Josephson effect, and the fact that high precision measurements are possible for such systems 
probably accounts for some of this interest. Another reason is that decay rates can be calculated 
theoretically based on models for such systems. The purpose here is to begin to introduce the 
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treatment of time dependence into these applications, and to discuss some of their quantum con- 
sequences. A perturbative time-dependent wavefunction will be developed for the time-dependent 
Schrodinger equation, and some applications will be considered. In particular, this type of wave 
function is applicable to the study of the Berry phase. 

The time-dependent Schrodinger equation in three space dimensions has the form 

ih^ = (1) 

where H is the Hamiltonian for the system. The case in which the Hamiltonian contains a potential 
function will be of interest here, in which case (1) can be written 

Here A is the Laplacian operator in terms of the position coordinates. To motivate the approach, 
consider a function %j} which is to satisfy (2) which has the particular form 

V>(x,t) = A(x,t)ei s ^. (3) 

In (3), both A and S depend on both the spatial variables and the time variable. Substituting (3) 
into (2), after differentiating with respect to the given variables, (2) takes the form 

' h m ~ A m = ~^ AA + 2 l VA ' vs + l AAS ~ /^ (VS)21 + VA - (4) 

Equating the real and imaginary parts on both sides of (4), two independent, coupled equations 
in terms of the variables A and S then result and are given by 

f + h^S? + V - = 0, (5) 

Ft A A 1 

^ + — AS + -VS- VA = 0. (6) 

ot 2m m 

Of course, by neglecting the term proportional to ti 2 in (5), the classical Hamilton- Jacobi equation 
in terms of the single function S results. In this approximation, a solution to (5) in (6) gives an 
equation for the function A. 

Let us now consider another version of this process. In other words, let if) have the basic 
structure given in (3), and then require that S satisfy the classical Hamilton- Jacobi equation. In 



fact, it proves more appropriate to introduce a function which is a superposition of such functions 
in powers of H. 

Consider a classical system which is governed by a Hamiltonian H, then the associated classical 
Hamilton- Jacobi equation has the following form in terms of r-space variables 

dS TTI dS dS . 

Although it should be possible to generalize what follows to higher dimensions, let us suppose 
that there is just one spatial variable in order to simplify the discussion. Then for a Hamiltonian 
of the form that generates (2), equation (7) takes the form 

2S<*> , + v + f = a (8) 

The Schrodinger equation (2) in one variable reduces to the form 

(*lfl2- „ + = (9) 



Define the linear operator 



L =^- v+i 4- < io > 



Equation (9) then takes the compact form 

L^ = 0. (11) 

Now the wave function ip in (3) will be generalized to the following series form in ascending powers 
of h by writing 

N 

*(x,t) = ^(^) fc a fc (x,t) exp(-^OM)). (12) 

fc=0 

In expression (12), the functions S(x,t) and the coefficients a,k(x,t) depend on both the space 
and time variables. Moreover, the a k are to be real-valued a k : K — > K. It will be seen that it is 
possible to force (12) to satisfy (11) up to a single error term which is proportional to a power of 
h by imposing a set of differential constraint equations on the coefficients a k which appear in (12). 
The second derivative of ^ is required in (9) and it is given by 



h 2 



I N N N 

^xx = ^—{-^^{ihfak^i^ + 2ihS x ^2(ih) k a ktX - (S x ) 2 ^2(ih) k a k 



2m 2m 

k=l k=0 k=0 



N s 

-(ih) N+2 a N;XX + ihS xx ^2(ih) k a k } exp(i— ). (13) 

k=0 

The time derivative is given by 

N „ 

im t = Y j m k+l ^M^) - S t *. (14) 

k=0 

Here, the variable subscript indicates differentiation with respect to x or t. Substituting these 
derivatives into Schrodinger equation (9), we obtain 

N 

2in ' ' ' 2/// ' ' /) 



h -* xx + iM? t = ^ihexp(^S) ^{-{ih) k a k _ hxx + 2S x (ih) k a Kx + S xx (ih) k a k } 

k=i 



1 i 1 i 

- — (ih) N+2 a N>xx exp(-S) + —ih{2S x a , x + S xx a } exp(-S') (15) 

1 ^ 
-(^* + V * + + £(^ fe+ V'exp(^). 

fc=0 

Now suppose that the function S(x,t) satisfies the Hamilton- Jacobi equation (8), then (15) takes 
the form 

1 i N 

L* = — exp(-S){^2(ih) k+1 (-a k _ hxx + 2S x a k , x + S xx a k + 2ma k>t )} 

k=l 

1 i 1 i 

+—ihexp(-S){2S x a , x + S xx a + 2ma , t } - — {ifr) N+2 a N , xx exp(-S). (16) 

It is now required that for a given function S, the coefficient functions a k must satisfy the following 
system of first-order partial differential equations with the definition a_i = 0, 

2S x ao x + S xx a + 2ma ,t = 0, (17) 

2^a fc)I + S xx a k + 2ma Kt - a k _ 1;XX = 0, k=l,---,N. (18) 

Equation (18) is non-homogeneous since the coefficient a k -\ and hence a k -\ jXX has been determined 
from (17)-(18) at the previous order k — 1. In the event that (17) and (18) hold, all the terms in 
the brackets on the right-hand side of (16) vanish, and the Schrodinger equation reduces to simply 

L * = -2^(^) jV+2 ^- ex P(^)" ( 19 ) 
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If wavefunction (12) had been an exact solution of (9), the right-hand side of (19) would have 
reduced to zero. Thus this procedure in which \1/ is taken in the form (12) provides a solution 
of (9) up to an error term proportional to h N+2 , provided that S is a solution of the classical 
Hamilton- Jacobi equation (8), and the a k satisfy the set of partial differential equations (17)-(18). 
Upon letting N — > oo, an asymptotic solution to equation (11) results which has the form 

oo 

^^(J2^(x,t)h k )eM^). (20) 

k=0 

Once S has been determined, the coefficients a k are obtained recursively as solutions of equations 
(17)-(18) which generate a time dependent solution. Convergence of (20) will not be discussed 
here. 

This procedure can be generalized to higher spatial dimensions by using the linearity of the 
Schrodinger equation (2). The Hamilton- Jacobi equation has the form 

d t S+^-(VS) 2 + V = 0. (21) 
2m 

Thus, if the wavefunction \l/ depends on three spacial coordinates x, the analogue of (12) will be 
written 

N 

*(x,i) = ^(^) fe a fc (x,t)exp(^(x,t)). (22) 

k=0 

The A\l/ term in (2) will include three expressions of the form (13) for each of the variables. 
Therefore, we can write 

h 2 

2m 

1 ' N 

= — ihexp(^S) J2l-( ih ) kAa k-i + 2(ih) k (VS) ■ Va k + (ih) k a k AS] 

k=l 

-(ih) N+2 Aa n exp(^) + -?-n[2(V5) • (Va ) + AS a ] exp(^) (23) 
n 2m n 

1 N 
~ { 2m~ {VS)2 + V + dtS) ^ + E^) fe+la ^ ex P(^)- 

k=0 

If we suppose that S is a function which satisfies the classical Hamilton- Jacobi equation (21), then 
the coefficient functions a k are required to satisfy the following system analogous to the pair (17) 
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and (18) 

2VS • Va + ASa + 2ma , t = 0, 

(24) 

2VS' • Va fc + ASa k + 2ma k , t - Aa k _ x = 0, k = 1, • • • , N 
Under the constraints (21) and (24), (23) reduces to 

m = --L(ih) N+2 Aa N exp(^). (25) 

Zillb lb 

Consider as an example the realistic case of the one-dimensional harmonic oscillator Hamilto- 
nian given by 

H = —p 2 + x 2 . (26) 
2m 

Then it is easy to check that the function 

S= J y/2m(P - x 2 ) dx - /3t (27) 

is a solution to (8) under the harmonic potential function. In (27), (3 is a separation constant 
obtained upon integrating (8). Substituting (27) into (17), a will be given by the first order 
equation 

2v / 2^v / ^ 3 ^$ £ " v 7 ^ . X a + 2m^ = 0. (28) 
dx y/p-x 2 dt 

This can be solved by the method of characteristics and the general solution is given by 

a (x, t) = (P- x 2 Y^(t - ^ arctan(^=L=)). (29) 

z \J p — x l 

Using (27), a in (28) can be written in terms of the derivative of S as 

ao(x,t) = (-=) 1 <p(t — arctan(— — ). 

V2m ^ o x 

In (29), <f is selected to be an arbitrary twice differentiable function of a single variable. The 

method of characteristics does not specify how ip is to be chosen, however, this can be done on 

physical grounds. Higher order coefficients for the expansion (12) can be obtained from (18) using 

(29), however they are rather more complicated and will not be presented here. 

As a final application of wavefunction (22), let us consider a curve C on a manifold of external 

parameters M, and the adiabatic evolution of the quantum system described by the parameter 



dependent Hamiltonian H = H(v) along the curve C, and v represents the set of parameters in 
vector form. A one-form may be defined on M as follows 



The integral of the one-form A^> along a curve C between v and v x in parameter space produces 
the following geometric quantity 



This is the Berry phase, where £ is an arbitrary two-dimensional submanifold of M such that 
<9£ = C and = dA<- n \ In terms of the wave function -0 n (x;v), the two-form is calculated 
from to be 



There is much interest in the study of the connection between the quantal Berry phase j n {C) 
and the classical Hannay angle. An important question in this regard asks if a classical system 
has a Hannay angle, will it also possess a Berry phase when it is quantized. Wavefunctions of 
the type discussed here can be useful for discussing questions regarding the relationship bewteen 
these quantities. 
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A in) = -Im(n\d M n). 
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